Control of scroll wave turbulence using resonant perturbations 
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Turbulence of scroll waves is a sort of spatio-temporal chaos that exists in three-dimensional 
excitable media. Cardiac tissue and the Belousov-Zhabotinsky reaction are examples of such media. 
In cardiac tissue, chaotic behaviour is believed to underlie fibrillation which, without intervention, 
precedes cardiac death. In this study we investigate suppression of the turbulence using stimulation 
of two different types, "modulation of excitability" and "extra transmembrane current". With 
cardiac defibrillation in mind, we used a single pulse as well as repetitive extra current with both 
constant and feedback controlled frequency. We show that turbulence can be terminated using either 
a resonant modulation of excitability or a resonant extra current. The turbulence is terminated 
with much higher probability using a resonant frequency perturbation than a non-resonant one. 
Suppression of the turbulence using a resonant frequency is up to fifty times faster than using a 
non-resonant frequency, in both the modulation of excitability and the extra current modes. We also 
demonstrate that resonant perturbation requires strength one order of magnitude lower than that 
of a single pulse, which is currently used in clinical practice to terminate cardiac fibrillation. Our 
results provide a robust method of controlling complex chaotic spatio-temporal processes. Resonant 
drift of spiral waves has been studied extensively in two dimensions, however, these results show for 
the first time that it also works in three dimensions, despite the complex nature of the scroll wave 
turbulence. 
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I. INTRODUCTION 

Turbulence of scroll waves is a sort of spatiotemporal 
chaos that is observed in some three-dimensional (3D) 
excitable media [l|, H, d, 0, [H, 0, 0, H, [§] ■ In cardiac tissue, 
such chaotic behaviour is known as fibrillation [To| and 
implies cardiac failure. 

The current method for terminating fibrillation in car- 
diac tissue is by means of a single electric pulse with a 
large amplitude. However, this approach is far from ideal. 
There are several known side-effects linked to the ad- 
ministration of the large electric shocks to patients [Til ] . 
Termination of fibrillation using shocks with a lower am- 
plitude would overcome such problems. 

Resonant drift of a spiral wave in two dimensions (2D) 
has been observed when one of the parameters of a model 
of an excitable medium was changed in time, with the 
period equal to that of the spiral wave [H, [r|. This 
phenomenon was later shown to be generic for reaction- 
diffusion excitable systems [HI, [HI- Thus it appeared 
that resonant drift could be used for moving the spiral 
wave around the medium to a boundary where it would 
terminate. 

Numerical experiments with reaction-diffusion models 
revealed that when close to boundaries, the period of 
the spiral wave changes, thus destroying the resonance in 
such a way that the drift trajectory turns away from the 
boundary. The main reason for this unruly behaviour is 
"resonant repulsion" , caused by the untuning of the res- 
onance between the spiral and the perturbation, due to 
the variation of the spiral's own frequency [3Q1I. This 
can be rectified by adjusting the frequency of the exter- 
nal forcing accordingly, based on some kind of feedback 



obtained from the re-entry itself. Feedback control of the 
resonant drift has been shown to overcom e rep ulsion from 
the boundaries and inhomogeneities [IH, Il6| . The same 
method can also eliminate multiple spiral waves, thus 
demonstrating that multiplicity of re-entrant sources in 
fibrillation is not in itself an obstacle for low-voltage de- 
fibrillation by this method [l5[ . 

Chambers of a heart, particularly ventricles, are 3D, so 
termination of scroll waves should be studied. Scroll wave 
turbulence presents a new challenge for resonant drift 
control, since here we are dealing not only with multiple 
sources, but sources that tend to multiply. 

A scroll wave rotates around a central filament. De- 
pending on the parameters of the medium, a filament of 
a circular shape may gradually contract or expand with 
time P, [l7|. For filaments of arbitrary shape, this prop- 
erty translates into "filament tension" Q . 

In excitable media, if a circular filament contracts, then 
a filament of any shape in the same medium has "pos- 
itive tension" and will shorten with time. Scroll waves 
with positive filament tension therefore either collapse or 
stabilise to a straight shape. In a bounded medium this 
can lead to the self-termination of the scroll. 

If a circular filament expands, then any filament shape 
is unstable as the filament has "negative tension" and will 
tend to lengthen. It was therefore conjectured P, 0, H| 
and subsequently demonstrated 0, 0, M, @ that such ex- 
citable media of sufficiently big size should support "tur- 
bulence" of scroll waves, where the scroll filaments grow, 
spontaneously bend, and break up to fragments upon 
collision with boundaries and with each other. Nega- 
tive filament tension is not the only mechanism of scroll 
wave turbulence: similar behaviour may occur due to 



non-uniform anisotropy of diffusivity, like that found in 
a ventricular wall[5j, [6| . 

Alonso et al. [3] considered the effect of applying a 
periodic non-resonant forcing on scroll wave turbulence 
produced by negative filament tension. By numerical 
simulations of the Barkley model, they showed that pe- 
riodic modulation of the medium's excitability with con- 
stant frequency higher than the frequency of the scroll 
waves can control the turbulence in the medium. They 
went on to propose a theory of this effect, based on the 
"kinematic description" of the scroll waves [18]. Their in- 
terpretation is that faster-than-resonant stimulation can 
effectively change the filament tension from negative to 
positive, thus disrupting the mechanism supporting the 
multiplication of scroll flaments, and compelling them to 
collapse. 

However, to our knowledge, the possibility of elimi- 
nating scroll turbulence by resonant stimulation has not 
been been investigated so far. We set this task for the 
present study. 

In this study we compare suppression of the scroll 
wave turbulence using (1) modulation of the medium's 
excitability, as used by Alonso et al. in 0], to enable 
comparison of our resonant forcing results with their non- 
resonant forcing, and (2) an "extra transmembrane cur- 
rent" forcing, as in [151 ]. Keeping in mind single shock 
cardiac defibrillation used in clinical practice, we also 
compare repetitive external forcing of constant and feed- 
back controlled frequencies to a single pulse extra current 
forcing. 

Our results show that resonant perturbation ensures 
the quickest termination of the scroll wave turbulence. 
Feedback controlled external perturbation was as effec- 
tive as constant frequency resonant perturbation, but of- 
fers the advantage of not having to know the correct fre- 
quency a priori. The resonant and feedback-controlled 
forcing suppress the 3D turbulence using amplitudes one 
order of magnitude lower than that of a single pulse cur- 
rently used in clinical practice to terminate fibrillation. 



II. METHODOLOGY 
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FIG. 1: (color online) Development of scroll wave turbulence 
from an initial scroll ring. The white lines show the filaments 
of the scroll waves. 



B. Numerical Methods 

For numerical simulations, we used EZscroll software 
by Barkley et al. [l|| [2(| , modified appropriately to de- 
scribe the stimulation. We used a 19-points finite dif- 
ference approximation of the Laplacian with equal dis- 
cretization steps Aa; in all three spatial directions, and an 
implicit first order Euler time stepping with a time step 
At. Simulations were run in a box (x, y, z) £ [0, L] 3 , with 
Neumann boundary conditions. In most simulations, we 
used Ax = 2/3, At = 1/30 and L = 60. 

In all simulations the model parameters were chosen 
as a = 1.1, b = 0.19 (or the average value of b(t) when it 
varied), and e = 0.02, as in Q. At this set of parameters 
a scroll wave will have negative filament tension. 

The choice of parameters was the same as in Q to 
allow comparison, with the exception of the discretiza- 
tion steps. We used cruder discretization steps, which 
allowed us to perform more simulations within reason- 
able CPU time. We performed also selected control sim- 
ulations with Ax = 0.4, At = 0.01 and L = 60, as in 0]; 
the results were quantitatively somewhat different but 
qualitatively similar (see below for details). 



A. Governing Equations 



The 3D numerical simulations presented here were per- 
formed using the Barkley model of excitable media [l9(, 
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where e is a small parameter e <C 1 characterising mutual 
time scales of the fast u and slow v variables, and a and 
b specify the kinetic properties of the system. Parameter 
b determines the excitation threshold and thus controls 
the excitability of the medium. The term h(t) represents 
an "extra transmembrane current" . 



C. Generation of turbulence 

The development of scroll wave turbulence is presented 
in fig. [TJ Starting at t = with the standard EZScroll 
scroll ring initial conditions in an unperturbed medium, 
the negative tension of the initial scroll ring caused elon- 
gation and bending of the filament. Interaction with the 
boundaries caused the filament to fragment and soon a 
complex tangle of many filaments filling the volume was 
observed. 

The turbulent state of the system was saved at five 
different times t = 240, 245, 250, 255, 260, and then each 
state used as an initial condition in our simulations. 
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D. Resonant Frequency 

There are different ways to define resonance between 
the forcing and the turbulence it is aimed to control. 
We considered three different frequencies, to which the 
forcing frequency can be compared: 

• The rotation frequency of an unforced vortex 
around its filament oj . It is also the frequency of a 
2D spiral wave in a large enough medium. It is typ- 
ically used as the leading-order approximation in 
any perturbative theoretical approaches, as 0,0- 

• The mean frequency of the unforced turbulence iZ>o . 
It is different from too due to interaction of scrolls 
with each other and with boundaries. This differ- 
ence can be significant as this interaction is the only 
factor that stops the filaments' growth in length 
and number. 

• The mean frequency of the forced turbulence 
Coq(A, cjf), which depends on the forcing amplitude 
A, and the forcing frequency ujj (sec fig. [2ja,b)). 
By definition, Hjq(0, ujf) = loq for any Uf. The dif- 
ference between uq and ujq is less obvious from the 
theoretical viewpoint than the difference between 
ujq and wo, as the theory of resonant drift of the 
scroll wave turbulence is yet to be developed. Yet 
we suppose that it is ujq that is to be compared to 
the forcing frequency to determine resonance, since 
it represents the de facto state of the controlled 
system regardless of the detailed mechanisms that 
brought it into that state. 

The frequency loq was measured for a single straight 
scroll. The frequencies (Do and u>o(A,u>f) were both mea- 
sured by recording the intervals Tj between the moments 
tj in which wavefronts passed through a recording point 
(x r ,y r ,z r ), that is, Tj = tj_^_i — t r , u(x r ,y r , z r ,t T j) = u*, 
where u r = 0, and ^f(x r ,y r , z r ,tj) > 0. The mean fre- 
quency for the entire simulation was then calculated as 
the average, <D = 2n/(Tj) = 2nN/Y,f=iTj, for all N 
intervals recorded, and similarly for u>q. For luq the sim- 
ulation was run for t € [0,5000]. For u)q(A, ojf) the sim- 
ulation was run for t £ [0, 5000] or until all scrolls were 
terminated if it happened sooner. 

For our crude discretization steps Ax = 2/3, At = 
1/30 and L = 40, when modulating the medium's ex- 
citability we have observed w = 0.90, <Do = 1.05 
and (Do (0.03, wo) = 0.74. The finer discretization steps 
Ax = 0.40, At = 0.01 and L = 60 produced w = 1-20, 
l)q = 1.27 and d)o(0.03, u)q) = 1.01. That is, the cruder 
discretiazation slows down the scroll waves overall, com- 
pared to the finer discretization, but the relationship be- 
tween the key frequencies remains similar. 



E. Forcing 

We investigated the application of the following types 
of forcing on scroll wave turbulence; 

1. Modulation of the medium's excitability, i.e. varia- 
tion of parameter b(t) around its average value bo- 

(a) Repetitive stimulation 
i. Constant frequency 

2. Extra transmembrane current, h(t) 

(a) Single pulse 

(b) Repetitive stimulation 

i. Constant frequency 

ii. Feedback controlled 

1. Modulation of the medium's excitability 

The medium's excitability in the model is defined by 
parameter b. Following 0, HH, we introduced into the 
model a spatially-uniform forcing by applying a periodic 
modulation of the parameter b in time, while keeping the 
extra current term zero, 

b(t) = b + Acos(oj f t), h(t)=0, (3) 

where bo = 0.19, A is the forcing amplitude, and u)f is 
the forcing frequency. Since b determines the excitation 
threshold, varying its value will vary the excitability of 
the medium. 

Starting from the five different turbulent initial condi- 
tions, the simulations were performed at different values 
of A and u>f, and the time taken for elimination of the 
turbulence was recorded. If any scrolls remained after 
t — 5000 then the experiment was stopped and consid- 
ered to have failed to eliminate the turbulence. 



2. Extra transmembrane current 

Here, simulations were performed with 

b{t) = b Q = 0.19, h(t) ± 0. 

(a) Single pulse. Simulations were started from the 
saved initial conditions at t = to and a single pulse of 
time duration A = 0.3 was applied, 

h(t) = AQ(t-t )Q(t + A-t), 

where 0() is the Hcaviside step function. After the shock, 
evolution of the filaments was observed for a further 250 
units of time. If no filaments remained at the end of this 
period, the shock was considered to be a success, and 
otherwise it was deemed to be unsuccessful. 

Shocks of different amplitudes A were tested, and the 
success threshold was defined as the amplitude which 
gives a 50% success rate over the five initial conditions 
used. 
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FIG. 2: (color online) The mean frequency £>o(A, u)f) of the perturbed turbulence, measured at the point (x r , y r , z r ) — (0, 0, 0), 
against forcing frequency u/ at different forcing amplitudes A due to (a) modulation of the medium's excitability and (b) 
extra transmembrane current forcing. In (a) and (b) the line wo = 1.05 is the mean frequency of the unperturbed turbulence, 
and the line cjo = 0.9 is the frequency of the single vortex. The resonant windows can be seen from the deviation of the 
mean frequency u>o(A,u>f) from the forcing frequency u)f due to (c) modulation of excitability and (d) extra current forcing. 
In (c) and (d) the line ojj — Qo(A, uif) — is drawn to highlight the resonant windows. The vertical brackets illustrate the 
corresponding resonant windows at forcing amplitude A = 0.03 for (c) and A = 0.3 for (d). 



(b) Repetitive stimulation. The stimulus h(t) was set 
to be a repetitive series of rectangular pulses of amplitude 
A, 

N 

h(t) = A^T0{t- t°) &(t s j + A - t). 

3=1 

We studied the effect of repetitive stimulation with both 
constant and feedback-controlled frequencies. 

i. Constant frequency 

Simulations were started from the saved initial con- 
ditions at t = ioj an d periodic pulses were applied, 
= t + j 2ir/ujf, j = 0,1,2.... Experiments 
were repeated for different forcing frequencies ujf 
and amplitudes A, and the time taken to eliminate 
the turbulence in each experiment was recorded. If 



any scrolls remained after t = 5000 then the exper- 
iment was stopped and considered to have failed to 
eliminate the turbulence. 

ii. Feedback controlled 

Stimulation h(t) was set to be a repetitive scries of 
rectangular pulses with timings t s - determined by 
taking feedback from the turbulence itself. Feed- 
back is taken from a recording point (x r ,y r , z r ) 
so that a pulse is applied every time that a wave 
passes through this point, that is, tj = t'j + ideiay, 
u(x r ,y r , z r ,t r j) =u* and ^{x r , y r , z r , t'j) > 0. We 
also varied the time delay ideiay for applying a pulse 
after a wavefront has passed through the record- 
ing point. Delays ranging from idciay = 0.0 to 
idciay = 6.0 were used, with increments of 0.1. 

Different locations were used for the recording point 
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FIG. 3: (color online) Different locations for the recording 
point 

(see fig. [3]): 

— in the corner of the domain, (x r , y r , z r ) = 
(0,0,0), 

— in the center of the domain, (x r ,y r , z r ) = 
(L/2,L/2,L/2), 

— in the center of a face of the domain, 
(x r ,y r ,z r ) = (i/2,L/2,0), 

— in the center of an edge of the domain, 
(x r ,y r ,z r ) = (L/2,0,0). 

Experiments were repeated for different locations 
for the recording point, different values of time de- 
lay £deiay, and amplitude A and the time taken for 
the elimination of the turbulence recorded. If any 
scrolls remained after t = 5000 time steps, then 
the experiment was stopped and considered to have 
failed to eliminate the turbulence. 



III. RESULTS 
A. Elimination of the turbulence 

Fig. [¥] illustrates the main result of our study: a res- 
onant stimulation can eliminate scroll wave turbulence, 
and docs it quicker and more reliably than a non-resonant 
stimulation. The figure shows evolution of the turbulence 
due to modulation of the medium's excitability at three 
different forcing frequencies at a fixed amplitude. The 
forcing frequency ujf = 0.8 is within the "resonant win- 
dow" (see below for a formal definition), and the turbu- 
lence is terminated quickly within t = 41. The forcing 
frequency u>f — 1.22 is above-resonant, and although the 
turbulence is terminated, it takes much longer, t = 2089. 
The forcing frequency u>f = 1.13 is also above the reso- 
nant window, and it leads to stabilisation of vortices in 
the center of the medium rather than their termination. 

We note here that the mechanism of "taming" of scroll 
wave turbulence suggested by Alonso et al. i s based 

on inversion of the filament tension from negative to pos- 
itive. The sequence shown in fig.[3]Jc) illustrates why this 



is not sufficient for defibrillation: the scroll filaments sta- 
bilise with a straight shape, which is consistent with their 
effective tension being positive, but it does not lead to 
their elimination. 

In the following subsections, we analyse in more details 
the empirical conditions required for successful termina- 
tion of the turbulence. 



B. Windows of resonant frequencies 

For both the modulation of the medium's excitability 
and extra transmembrane current forcing, we have varied 
the frequency u) / and amplitude A to assess their effects 
on termination of the turbulence. We observed very dif- 
ferent effects of the modulation of excitability and of the 
extra current forcing on the mean frequency of the turbu- 
lence LJo(A,ujf). In this section we define a resonant win- 
dow of frequencies for each amplitude and forcing type. 

1. Modulation of excitability 

As the amplitude of the modulation increases, the 
mean frequency of the turbulence u>o(A,u)f) decreases 
dramatically (see fig. [2ja)). When using the largest forc- 
ing amplitude (A = 0.05), the frequency of the turbu- 
lence reduced to over half that of the frequency of the 
unperturbed turbulence. 

Resonant windows can be identified in fig. H^c). We 
define the resonant window to be the range of forcing 
frequencies ujf, for which uif w u)o(A,u)f). The upper 
and lower bounds for the resonant window can be seen 
as the first and last points where u>f — uio(A, uif) = 0. The 
resonant window is taken to be this range and a further 
0.01 either side of this range. 

There is a different resonant window for each ampli- 
tude. As the amplitude increases, the size of the reso- 
nant window increases and shifts towards lower forcing 
frequencies. 

The above definition of the resonant window should 
be used with caution for the lowest amplitudes. E.g. for 
A = 0.01 in FigJ^c) the window must be between 0.873 
and 1.045. However, most of this interval corresponds 
to a "false resonance" , when uif ks ujq but that does not 
lead to termination. Termination of the turbulence at 
this forcing amplitude is observed in the narrow vicinity 
of ujf = 0.873 only. Applying a forcing with frequency in 
the vicinity of uif = 1.045 maintains the turbulence. 

2. Extra transmembrane current 

As the amplitude of the extra current forcing increases, 
the mean frequency of the turbulence uio(A,uif) decreases 
(see FigJSJb)), though not so dramatically as in the 
case of modulation of the medium's excitability, compare 
FigEI a ) an d (b). Even for the largest forcing amplitude 
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A = 0.5, the reduction in the frequency of the turbulence 
due to the extra current forcing is not more than 20% of 
the frequency of the unperturbed turbulence. 

Resonant windows for extra current forcing can be 
identified in Fig|2jd). The resonant windows were de- 
fined in the same way as in section IIII B 11 There is a 
different resonant window for each amplitude, although 
for all forcing amplitudes that we tested, the resonant 
windows are in the vicinity of the frequency of a straight 
scroll wo = 0.9. 

For the lowest amplitude A = 0.1 in Fig[2^d) there 
is no obvious resonant window. However, the resonant 
termination of the turbulence at this forcing amplitude 
is observed in the narrow vicinity of 0.870. The interval 
of frequencies above that and up until 1.09 corresponds 
to the "false resonance" . 



C. Termination times 

1. Modulation of the medium's excitability 

Fig. [5] presents the turbulence termination times due 
to modulation of the medium's excitability at four differ- 
ent forcing amplitudes. The vertical brackets show the 
windows of resonant frequencies for each amplitude, as 
in fig. [21 From these data, it can be seen that termina- 
tion of turbulence is fastest when the forcing frequency 
ujf is within the resonant window. The turbulence can be 



terminated with a frequency outside of the resonant win- 
dow. Although, further away from the resonant window, 
average termination time increases and the probability of 
success decreases. 

Fig.[6]shows a plot of termination time against (a) fixed 
forcing amplitude A, and (b) untuning of the resonance 
S defined as: 

s = \uf -u) {A,u) f )\ 
a 

where a is the standard deviation of the turbulence fre- 
quency recorded throughout a simulation, 

The strength of resonance is a measure of how close the 
forcing frequency is to the frequency of the scroll wave 
turbulence. 

It can be seen from fig. Ufa) that increasing A reduces 
the termination time. The reduction is not very pro- 
nounced at larger A; one should bear in mind here that 
the data in this graph are for all frequencies, resonant 
or not. Fig. [BJb) shows that at smaller S, i.e. a better 
resonance, the time taken to eliminate the turbulence re- 
duces. Here the data are for all forcing amplitudes, large 
and small. Comparing fig. E^a) and fig.[6^b), and taking 
into account that reliable forcing amplitudes for the the 
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(c) A = 0.03 ^ (d) A = 0.05 * 

FIG. 5: (color online) Modulation of medium's excitability: termination times at different amplitudes: (a) A = 0.01, (b) 
A — 0.02, (c) A — 0.03 and (d) A = 0.05. Black dots: termination times for individual simulations. Red solid line: median 
values of the termination times at every fixed frequency. Blue dashed line: same, geometric mean values. For averaging and 
visualization purposes, we assign value t = 10 to the failures. Black vertical brackets: the windows of resonant frequencies for 
each amplitude. 




FIG. 6: (color online) Modulation of medium's excitability: Termination time against (a) fixed forcing amplitude A, (b) 
untuning of the resonance S. The dots represent termination times for individual simulations. The solid line goes through the 
median values of the termination times at either (a) or (b), and the dashed line through the geometric mean values. 
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FIG. 7: (color online) Single pulse: termination of the turbu- 
lence at forcing amplitude A = 4.3. 



turbulence termination seems to be A > 0.02, termina- 
tion times are more sensitive to the quality of resonance 
than to the forcing amplitude. 



2. Extra transmembrane current forcing 

(a) Single pulse. The single pulse stimulation was 
tested having in mind the current clinical practice for 
cardiac defibrillation is by means of a single electric pulse 
of large amplitude. We defined the single pulse success 
threshold as the amplitude at which the turbulence is 
terminated in more than 50% of the experiments. In our 
setup, this success threshold was found to be A = 4.3. 
Fig. [7] shows an example of successful defibrillation with 
a single pulse shock at amplitude A = 4.3. 

(b) Repetitive pulses. Fig. [H] presents the turbu- 
lence termination times, both for constant-frequency and 
feedback-controlled extra currrent forcing. The solid- line 
and dashcd-line curves show the dependence of the ter- 
mination times on the extra current forcing frequency u>f, 
at fixed amplitude A = 0.3. The vertical brackets show 
the window of resonant frequencies as in fig. Hfd). The 
dashed horizontal straight line shows the geometric mean 
termination time for the feedback controlled stimulation, 
and the solid horizontal line shows the corresponding me- 
dian termination time. From these data, it can be seen 
that forcing frequencies within the resonant window en- 
sure the fastest termination of turbulence. The turbu- 
lence can be terminated with a frequency outside of the 
resonant window but there the probability of the turbu- 
lence termination decreases. The individual experiments 
with failed termination are depicted as log 10 (i) = oo, 
and counted for the sake of averaging as log 10 (i) = 5. 
Further away from resonant frequencies the turbulence 
termination time rapidly increases. 

Fig. [5] illustrates the effects of the stimulation ampli- 
tude A. The left panels in it are similar to fig. [8] and 
present the turbulence termination time dependence on 
the extra current forcing frequency u>f, at four differ- 
ent amplitudes A = 0.1,0.2,0.4,0.5. The histograms on 



the right show the distributions of termination times for 
feedback controlled experiments, for the four different lo- 
cations of the recording point. All the observations made 
for fig. [5] are valid for the forcing amplitudes in fig.O 

For amplitudes A > 0.2 in feedback controlled ex- 
periments, the average termination time is close to the 
average termination time achieved at resonant frequen- 
cies. For lower amplitudes, the success probability using 
feedback-controlled controlled stimulation falls down, in 
the same way as it does for the constant-frequency forc- 
ing. 

The location of the recording point used for the feed- 
back is also important for a successful termination. The 
most successful locations, in the 3D experiments, appear 
to be the corner or the edge of the medium. The center 
location for the recording point appears to be the worst 
for all forcing amplitudes tested. 

Fig. [TO] shows evolution of the turbulence due to ex- 
tra current forcing at different forcing frequencies ojj = 
0.9, 1.22, 1.15, and the feedback controlled, at the fixed 
amplitude A = 0.3. The forcing frequency cj/ = 0.9 is 
within the resonant window, and the turbulence is ter- 
minated quickly by t = 134 (series a). The feedback con- 
trolled stimulation terminates the turbulence by t = 190 
(series b). The forcing frequency uif = 1.22 is above- 
resonant, although the turbulence is terminated it takes 
ten times longer, to t = 1773 (scries c). The forcing 
frequency u>f = 1.15 is also above the resonant window. 
Simulation with that frequency leads to stabilisation of a 
vortex in the center of the medium (series d) , similar to 
what was observed for faster-than-resonant modulation 
of excitability, sec fig-HIc). 

Fig- SH is similar to fig. [6] and shows plots of termina- 
tion time against (a) fixed forcing amplitude A, and (b) 
untuning of the resonance 5, defined in the same way as 
for the modulation of excitability forcing. It can be seen 
from fig. Illf a) that increasing the amplitutde A reduces 
the termination time. Fig. mT b) shows that reducing the 
untuning S also reduces the termination time. Compar- 
ing fig. mT a) and fig-E^b), and taking into account that 
reliable forcing amplitudes for the turbulence termination 
seems to be A > 0.2, termination times are more sensi- 
tive to the untuning of resonance (5 than to the change of 
forcing amplitude A. 



IV. DISCUSSION 

A. Resonant stimulation terminates scroll 
turbulence 

We have shown that a termination of the scroll wave 
turbulence can be achieved by repetitive stimulation. 
This is despite the fact that scroll waves continue to grow 
and multiply while we drive them to elimination. We 
have studied two different methods of forcing, the modu- 
lation of medium's excitability and the extra transmem- 
brane currents forcing. Where comparable, the qualita- 
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FIG. 8: (color online) Extra current forcing: termination time for the amplitude A = 0.3. The black dots represent termination 
times for individual simulations. The red solid line goes through the median values of the termination times at a fixed frequency, 
and the blue dashed line through the geometric mean values. The horizontal straight lines show the median (red solid line) 
and geometric mean (blue dashed line) termination times for the feedback experiments. The vertical brackets designate the 
window of resonant frequency. 



tive depenencies for the two methods were similar. For 
a successful termination, the amplitude of the repetitive 
forcing should be higher than a certain threshold. How- 
ever, this threshold is still much lower than that required 
for termination by a single shock. The termination is 
achieved with the highest probability, and in the quickest 
time, when using a resonant forcing frequency, i.e. when 
the frequency of stimulation, u>f, is close to the de facto 
frequency of the forced scroll waves, u>o(A,ujf). Namely, 
we have shown that for both types of forcing, the turbu- 
lence termination becomes faster for smaller values of the 
untuning of the resonance 5 (see fig.EJb) and fig. fTTT b)). 



B. Resonant windows 

We have found that the resonance between the forcing 
and the scrolls is characterised not by a single resonant 
frequency, but by a resonant window which depends on 
the type of forcing and its amplitude. For the modula- 
tion of the medium's excitability, larger amplitudes cor- 
responded to wider resonant windows, shifted towards 
lower forcing frequencies. For the extra transmembrane 
current forcing, the resonant windows did vary in size and 
location, but for the amplitudes that we tested they re- 
mained close to the frequency of a single straight vortex, 
ujq « 0.9. 



Our results also show that for both modulation of ex- 
citability and extra current forcing, the fastest termina- 
tion was achieved when the forcing was applied with a 
frequency chosen within the resonant window. Termina- 
tion using a forcing frequency outside of the resonant win- 
dow occurs with a much lower probability and a longer 
termination time. 

Existence of resonant windows rather than unique reso- 
nant frequencies may be a purely statistical phenomenon 
due to fluctuations of frequencies of the scroll turbu- 
lence, or may be an indication that the forced turbu- 
lence adjusts its frequency in response to the forcing, i.e. 
a "frequency locking" in average (perfect locking is not 
feasible between a periodic forcing and chaotic turbu- 
lence). Nothing like this has been reported for resonant 
drift in 2D, thus we may be dealing with a specifically 
three-dimensional phenomenon. Indeed, 3D scrolls have 
additional degrees of freedom compared to 2D spirals, 
e.g. twisted scrolls can rotate faster than straight scrolls 
[2ll . [22| ■ This might offer an explanation of cases of false 
resonance, when formally defined resonant windows are 
abnormally extended towards higher frequencies but are 
not associated with fast and reliable termination, see e.g. 
the case of A = 0.01 in fig. [UJa,c) and fig. [5fa). Indeed, 
in case of strong twist, different parts of the same fil- 
ament have different phases and are forced in different 
directions, which does not result in an overall directed 




FIG. 9: (color online) Extra current forcing: termination times for the amplitudes (a) A — 0.1, (b) A — 0.2, (c) A — 0.4, 
(d) A = 0.5. The left panels are as in Fig(8] The histograms on the right are distribution of termination times for feedback 
controlled experiments, at different locations for the recording point (they are rotated 90° clockwise and flipped in the vertical 
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FIG. 11: (color online) Extra current forcing: Termination time against (a) fixed forcing amplitude A, (b) untuning of the 
resonance S. The dots represent termination times for individual simulations. The solid line goes through the median values of 
the termination times at either (a) or (b), and the dashed line through the geometric mean values. 



(a) 



A m pi it u dc 



(b) 



12 



movement and does not bring about termination. 



C. Resonant vs non- resonant stimulation 

It has been previously shown Q in Barkley's model 
with the same model parameters as we used here, that 
scroll wave turbulence can be controlled by a weak non- 
resonant modulation of the medium's excitability. A the- 
ory was presented in [l8| explaining that this control of 
turbulence was due to an inversion of the filament ten- 
sion from negative to positive, which can happen if the 
frequency of forcing is higher than the frequency of the 
scrolls. It was argued that such stimulation causes the 
filaments to collapse and could therefore be used for ter- 
mination of the scroll wave turbulence. We have seen, 
however, that positive tension may lead to stabilization 
rather than termination of scrolls, sec fig-BJc) and lTUT d). 

In Q an above-resonant frequency forcing was used 
to control the turbulence. More specifically, their forc- 
ing frequency u>f was almost equal to (1% higher than) 
the frequency of a straight scroll, u>q. The frequency of a 
forced turbulence is significantly lower than the frequency 
of a straight scroll, e.g. u>o{A, loq) < u>o, thus forcing with 
frequency loq is above-resonant. This is true both for the 
finer discretization steps used in Q and cruder discrteti- 
zation steps used in a majority of our simulations. 

To make a specific comparison, let us consider ex- 
citability modulation with amplitude A = 0.03 (same as 
in Q) and frequency uif = 0.91 which is about 1% higher 
than uiq. As can be seen from fig. 0(c) , such forcing gives 
a mean turbulence termination time of t = 251 (compare 
with termination time t = 1510 in the example shown 
in Q). Within the resonant window for this amplitude 
as defined in this article, the mean termination time is 
between t = 71 and t = 32, i.e. 3 to 7 times faster than 
using above-resonant forcing frequency as in 0]. 

So, a direct like-for-like comparison shows that al- 
though the above-resonant frequency stimulation sug- 
gested in 0] works in principle, the resonant stimulation 
works more reliably and much faster. 

It has been reported by Wu et al. [23| that using a 
travelling-wave modulation of the mediums excitability 
can control scroll wave turbulence faster than the modu- 
lation used in Fig. 6 in [23| clearly shows that the op- 
timum forcing frequency is below the rotation frequency 
of an unforced vortex uq. Alas, Wu et al. did not con- 
trol the de facto frequencies of the scroll turbulence so it 
is not possible to interpret their results unambiguously. 
However, as we have shown here that the resonant fre- 
quencies are below loq, it is quite possible that the real 
reason for the advantage achieved in [23|] compared to Q 
is not (only) in using travelling waves, but simply in us- 
ing a frequency within a resonant window. A definitive 
answer to this question requires further investigation. 



D. Feedback control works in 3D 

As the resonant window may not be known a priori, we 
tested a simple algorithm for the feedback control with 
the extra current forcing. It has been previously shown in 
2D that applying a repetitive feedback controlled forcing 
causes a spiral wave to drift to a boundary along a pre- 
dictable trajectory [lH. fl^. and possibly terminating it 
faster than with the constant frequency forcing, as feed- 
back forcing can overcome the "resonant repulsion" of 
the drifting spirals from boundaries and inhomogeneities 
fl4| . Simulations of spiral waves in the two-dimensional 
version of our model easily demonstrate resonant repul- 
sion (see fig. [T2T a)). so it must play some role in 3D 
behaviour as well, even though it may be not straightfor- 
ward in scroll turbulence. Our 3D experiments show that 
the termination with a feedback-controlled forcing is nev- 
ertheless achievable even at relatively weak amplitudes. 

Our results also show that the location of the record- 
ing point is important for a successful feedback-controlled 
termination, which is in good agreement with earlier ob- 
servations of feedback-driven resonant drift in two dimen- 
sions [U Hfl. The most successful locations, in the 3D 
experiments, appear to be the corner or the edge of the 
medium. The center location for the recording point ap- 
pears to be the worst for all forcing amplitudes tested. 
It has been shown in 2D experiments [26| that a line of 
recording points is a robust approach. Therefore, if the 
same holds for 3D, a line of recording points down one 
edge of the medium may be the optimal choice. 



E. Feedback controlled vs constant frequency 
stimulation 

Typically, the average termination time using feedback 
is close to that achieved at resonant constant frequencies. 
At the lowest stimulation amplitudes, both feedback- 
controlled and constant-frequency termination times in- 
crease, and the probability of success decreases. 

For a range of frequencies within or near the reso- 
nant window, the termination was achieved, on aver- 
age, quicker than with feedback. To elucidate a pos- 
sible reason for this difference, we have considered a 
two-dimensional version of the Barkley model, with the 
same model parameters and applied the same feedback 
controlled forcing as in our 3D simulations. We have 
found that applying a forcing with feedback in this case 
causes the spiral wave to drift along a rather compli- 
cated "snaky" trajectory, shown in fig. [T27 b). This be- 
haviour is similar to that described recently by Zykov et 
al. [13, HI] , and is caused by an instability, related to the 
delay between a change of the position or phase of the 
spiral and its detection by the feedback electrode, due 
to the distance between them and a finite speed of the 
waves. Indeed, in fig. [T2]the spiral core is more than one 
wavelength away from the recording point. It can be seen 
from fig. [T27 b) at t — 525 that this resonantly drifting 
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FIG. 12: (color online) Two-dimensional phenomena, (a) 'Resonant reflection': Extra current forcing with constant frequency 
ujf = 1.1855 and amplitude A = 0.01. (b-d) 'Snake instability': Extra current forcing with feedback controlled frequency 
applied with amplitudes A = 0.01 (b), A = 0.004 (c) and A = 0.001 (d). Model parameters a = 1.1, b = 0.19, e = 0.02, L = 60, 
Ax = 0.4 and At = 0.01, with 9-point approximation of the Laplacian. The recording point in (b-d) was located in the bottom 
left corner. 



spiral does not terminate at the boundary. Instead, it 
embarks on a continuous loop near the boundary of the 
medium. 

According to [U [28J, this 'snake instability' should 
disappear at lower amplitudes. Indeed, this is what hap- 
pened in our simulations, sec fig. I12f c.d). The instabil- 
ity was less pronounced when the amplitude is decreased 
from A = 0.01 in (b) to A = 0.004 in (c), and was com- 
pletely gone for A = 0.001 in (d). However, although 
the drift trajectory towards the boundary was shorter, 
the drift velocity was smaller proportionally to A, so the 
time to reaching the boundary did not decrease. Besides, 
the decrease in A made annihilation at the boundary even 
less likely: in (c) the spiral stuck near the boundary sim- 
ilar to (b), whereas at a further reduction of A in (d) it 
embarked on a drift along the boundary which was faster 
than its resonant drift in the center of the domain. So, 
reducing the amplitude is not necessarily a satisfactory 



solution to the problem of the 'snake instability' interfer- 
ence with terminating the vortex. Alternative solution 
could be to reduce the delay in the feedback, say by us- 
ing a global (ECG) rather than local (point electrogram) 
signals. However, in known computer models and exper- 
imental observations, scroll wave are apparently rather 
large in the scale of the heart chambers and are typically 
no more than one wavelength away from a boundary. 
Thus the 'snake instability' may not be a problem in a 
real heart. 
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